Generation of multi-photon Fock states by bichromatic adiabatic passage: 
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We propose a robust scheme to generate multi-photon Fock states in an atom-maser-cavity system 
using adiabatic passage techniques and topological properties of the dressed eigenenergy surfaces. 
The mechanism is an exchange of photons from the maser field into the initially empty cavity by 
bichromatic adiabatic passage. The number of exchanged photons depends on the design of the 
adiabatic dynamics through and around the conical intersections of dressed eigenenergy surfaces. 
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I. INTRODUCTION 

Over the past few years, new sources of anti-bunched 
light that are able to emit a single photon in a given time 
interval has been the subject of an intense theoretical 
and experimental research. The driving force behind the 
development of these non-classical sources is a range of 
novel applications in quantum information theory which 
builds on the laws of quantum mechanics to transmit, 
store, and process information in varied and powerful 
ways. Advances in this field rely on the ability to manip- 
ulate coherently isolated quantum objects while eliminat- 
ing incoherent interactions with the surrounding environ- 
ment. Single photon states act as elementary quantum 
bits (qubits) in quantum cryptography Q, |j| and tele- 
portation of a quantum state |j] where their entangled 
states enable secure transmission of information. 

Many different types of single-photon sources have 
been proposed and realized using the controlled excita- 
tion of single molecules or of single nitrogen- vacancy 
centers in diamond nanocrystals 0, controlled injection 
of carriers into a mesoscopic quantum well |^ and using 
pulsed excitation of semiconductor quantum dots ■ 

In the context of cavity QED, single-photon Fock 
states have been produced by a Rabi 7r-pulse in a mi- 
crowave cavity [ll and by the STIRAP technique in 
an optical cavity I3 based on the scheme proposed in 
[T^ where the Stokes pulse is replaced by a mode of a 
high-Q cavity. The STIRAP process has also been stud- 
ied in a system of four-level atom interacting with a cav- 
ity mode and two laser pulses, with a coupling scheme 
which generates two degenerate dark states [l^. In all 
these cavity QED schemes one atom interacts with a 
single-mode high-Q cavity and generates one photon. As 
the atoms pass through the cavity one by one. more pho- 
tons can be added to the cavity. Recently another scheme 
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has been proposed to generate a two-photon Fock state 
by a single two-level atom interacting with a supercon- 
ducting cavity which sustains two non-degenerate orthog- 
onally polarized modes . The photons are transferred 
from the source mode into the target mode of the cavity 
by a third-order Raman process. However this scheme is 
not robust relative to the velocity of atoms. 

In this paper we propose a scheme in which a two- 
level atom interacts counterintuitively with a single- 
mode high-Q cavity and a delayed maser field that are 
both near-resonant with the atomic transition, allowing 
to produce a controlled number of photons in the cavity 
depending on the design of the adiabatic passage. This 
process is referred to as a bichromatic adiabatic passage 
since two near-resonant interacting fields act on a single 
transition. A related work involving exchange of pho- 
tons between two laser fields through a bichromatic pro- 
cess can be found in Ref. The transfer of photons 
from the maser field into the cavity field is based on the 
adiabatic passage between two dressed states which are 
the eigenstates of the coupled atom-maser-cavity system. 
This process is robust because it does not depend on 
the precise velocity of the atom or on the precise tuning 
of the maser and the cavity frequencies. The dynam- 
ics of the process, under the adiabatic conditions, can 
be described completely by the topology of the dressed 
eigenenergy surfaces. This topological aspect is the key 
to the robustness of the process. Our method is based 
on the calculation of the dressed eigenenergy surfaces of 
the effective Hamiltonian as a function of the two Rabi 
frequencies associated to the maser and the cavity fields, 
and the application of adiabatic principles to determine 
the dynamics of the process in view of the topology of 
the surfaces. 

The paper is structured as follows. In Sec. II, we 
use the Floquet formalism and the phase representation 
of the creation and annihilation operators to construct 
the effective Hamiltonian of the atom-maser-cavity sys- 
tem. Eigenenergy surfaces of the effective Hamiltonian 
are displayed in Sec. HI function of the normal- 
ized Rabi frequencies of the cavity and the maser fields. 
We demonstrate how the analysis of these surfaces al- 
lows to design different adapted adiabatic paths leading 
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FIG. 1: Experimental configuration and level scheme of the 
atom. 



to different photon transfers into the cavity field without 
changing the atomic population at the end of the inter- 
action. Sec. IV is devoted to the numerical simulation of 
the evolution governed by the effective Hamiltonian and 
the final probabilities of the one and two-photon transfer 
states. Finally, in Sec. V we give some conclusions and 
indicate conditions for experimental implementation. 



II. CONSTRUCTION OF THE EFFECTIVE 
HAMILTONIAN 

We consider a two-level atom of upper and lower states 
and |— ) and of energy difference E+ — E- = ojo as 
represented in Fig. 1. We use atomic units in which 
h = 1. The atom in its lower state is released from a 
source of atoms and falls through a high-Q cavity with 
velocity v. The atom first encounters the vacuum mode 
of the cavity with frequency ujc and waist Wc and then 
the maser beam with frequency wm and waist Wm- Both 
the maser and the cavity fields are near-resonant with 
the atomic transition. The distance between the crossing 
points of the cavity and the maser axis with the atomic 
trajectory is d. The travelling atom encounters time de- 
pendent and delayed Rabi frequencies of the cavity and 
the maser fields: 



G{t) = 

n{t) = 



UJC 



= -/if. 



(1) 



where /i, Vi„odc, are respectively the dipole moment of 
the atomic transition, the effective volume of the cavity 
mode and the amplitude of the maser field. The detuning 
of the maser and the cavity fields from the atomic transi- 
tion arc Ai = wq— Wi, i = C, M . We take the frequencies 
of the fields such that their difference 5 = ojc — is 
positive and very small with respect to ujc and lom . Also 
we assume 



max{|G(t)| , <tJo. 



(2) 



Under these conditions, the counter-rotating terms can 
be discarded in the rotating- wave approximation (RWA) . 



The semiclassical Hamiltonian of the atom-maser-cavity 
system can thus be written as 



H{t) = ujca^a I2 + 



LOq 



G{t) 



a 
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(3) 



where a, are the annihilation and creation operators 
of the cavity field, I2 is the 2x2 identity matrix and 
the phase Oq is the initial phase of the maser field. The 
energy of the lower atomic state has been taken as 0. 
This Hamiltonian acts on the Hilbert space TL®J- where 
7i = is the Hilbert space of the atom generated by 
±) and T is the Fock space of the cavity mode generated 
by the orthonormal basis {|n); n = 0,l,2,---} with n 
the photon number of the cavity field. The dynamics of 
system is determined by the Schrodingcr equation 



(4) 



where (j){t) G Ti. (E) J- with initial condition (/)o = (j){to) = 
|— ) (8) |n = 0). We can think of Eq.Q as a family of 
equations parameterized by the phase 60 . The time- 
dependent Hamiltonian ^ contains two different time 
scales: the period T = 2tt/uj characterizing fast oscilla- 
tions of the maser field and Ti„t w ~ character- 
izing the slow change of the field amplitudes G{t),^l{t). 
The fast periodic time-dependence can be taken into ac- 
count by use of Floquet theory. The evolution equation 
in the Floquet representation reads: 



(5) 



where 9 appears as a dynamical variable defined on the 
unit circle §^ of length 2tt, and K(t-, 9) is the Floquet 
Hamiltonian 



K{t;9) = ujca)a I2 



LUq 





Git) 
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n{t)/2 







,(6) 



The only time dependence of the Floc^uet Hamiltonian 
is from the slow variation of the field amplitudes. This 
Hamiltonian acts on the enlarged Hilbert space ^ = 
U®T ® H where H := L2(S\ d6'/27r) denotes the space 
of square integrable functions on the circle §^ of length 
27r, with a scalar product 



d9 
2^ 



riio)f2{0). 



(7) 



This space is generated by the orthonormal basis 
{e''^^ ; fc G Z}. One can interpret k as the relative pho- 
ton number with respect to the (large) average photon 
number k of the maser field. The operator —i-^ can 
be interpreted as the relative photon number operator 
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of the maser field [20j. The eigenvectors of the zero- 
field Floquet Hamiltonian are |±, n, k) = |±) ® |n) ® e*''^ 
which form an orthonormal basis of the enlarged Hilbcrt 
space /C. The relation between the solutions of Eqs. Q 
and Q is established as follows If ip{t]9) is a solu- 
tion of (0) with initial condition "0(^0; 9) = 4>o<E) Ic, then 
(j){t; 6q) := ip{t; Oo+uit) is a solution of with the initial 
condition 0(io) = 0o- l£ is the basis function e'*^^ with 
fc = 0. We remark that if at the end of interaction t — tf, 
the solution of lO has a form ip^tf^O) = 0/ (8) e^^^ then 
the probability for the solution of Q to be found in the 
final states |±,n), i.e. \{4){tf;0o)\±,n)\'^ = |(0/|±,n)p, 
will not depend on the phase OQ+wtf of the semiclassical 
Hamiltonian JSJ. 

The evolution of (0 due to slow field amplitudes will be 
treated in the enlarged Hilbcrt space by adiabatic prin- 
ciples. We first show that the dynamics of © under 
the bichromatic interaction can be described by an effec- 
tive Hamiltonian. We start by applying to the Floquet 
Hamiltonian the unitary transformation 



R = 



(8) 



which yields 
K' 



R^KR ^ LOca)a I2 - I2 

oO 

( Am nit)/2\ f ae^' 



(9) 



The third term of K' , denoted HjiwA, is the so-called 
RWA Hamiltonian, associated to the maser field and the 

atom. Its eigenvalues are 2Aj!^'' = Am ± \/ {Am)^ + {^)^- 
To simplify and decouple the Hamiltonian (|5J), we use 
the phase representation of a and as formulated by 
Bialynicki-Birula 



■ 9 t 

dip 



■ 9 +iv t ■ ^ 

e ^ , a'a ^ —i 



dip 



which gives 

, d d 



dip 
(10) 



G{t) 











Defining the new variables 

C := ^ - 0, ^l := 0, 

we have 

d d d d d 

dip dC 99 drj d( 

The eigenbasis of {-i-^ - is {e^'Ve^ke . ^ ^ 

0, 1, 2, • • • , = 0, ±1, ±2, • • • } which can be written as 



(12) 
(13) 



e 



e -e 



(14) 



where m := n + fc = 0, ±1, ±2, • ■ • . Substituting H13|l in 
((TT|l gives 

d d 



G[t) 



(15) 



We can define new operators as 




(16) 



which verify the standard commutation relations [6, 6^] = 
1. The new bosonic operator h that corresponds to the 
process of creation of a cavity photon and associated an- 
nihilation of a maser photon, can be intuitively inter- 
preted as the transformation of a maser photon into a 
cavity photon. The Hamiltonian 1)15(1 can thus be ex- 
pressed as 



K' =-iWM^t2+H-'', 

di] 



where H is the reduced effective Hamiltonian 



(17) 



( Am ^^(t)/2\ 
\^{t)/2 ) 



Git) 



b 
&t 



(18) 

K' is defined on the Hilbert space generated by the or- 
thonormal basis {|±)®e"'^(X)e™'' ; n = 0, 1, 2, • • • ; to = 
0, ±1, ±2, • • • } and H'^^ is defined on the Hilbert space 



generated by the orthonormal basis {|±) 
0, 1, 2, • • • } where n is the number of exchanged photons 
from the maser field into the cavity field. 



III. TOPOLOGY OF THE DRESSED 
EIGENENERGY SURFACES 

The dressed eigenenergy surfaces of K' TT^ can be 
calculated numerically and can be displayed as a function 
of the normalized Rabi frequencies G/6 and ^1/5. These 
surfaces are grouped in families for different values of 
TO, each of which for zero fields consists an infinite set 
of eigenvalues with equal spacing Am- In what follows, 
we study the to = family only which means k = —n 
(fc £ Z is the relative photon number of the maser field 
and n > is the photon number of the cavity field). 
The labelling of the dressed eigenenergy surfaces can be 
performed in terms of the eigenvectors of the zero-field 
original Hamiltonian, 



K'{n^O,G = 0) = {Sb^b 



drj 



Am 




(19) 
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with eigenvalues 



£;;,„^„(f^ = 0,G = 0) = 6n 
E'_„,^{n^Q,G = 0) = 6n 



mujM 
mujM . 



(20) 



Since the eigenvectors of K and K' are related by the 
transformation JHJ as = R\(p)' , the correspondence 
between the eigenvalues of the zero-field effective Hamil- 
tonian and the eigenvectors of the original zero-field 
Hamiltonian are 



E'in^o,G^o)+^„,m=o ^ \+,ny 

^'(r!-0,G-0)„,„,™=o ^ \-,ny 



+ , n, —n — 1) 
\-,n,-n). 

(21) 



Figure 2 represents the m = family of the eigencn- 
ergy surfaces of K' as function of the instantaneous nor- 
malized Rabi frequencies G/S and Q/S. Any two neigh- 
boring surfaces have conical intersections on the plane 
G = and also on the plane = (except the first sur- 
face), corresponding to the situations where only one of 
the fields (maser or cavity) is interacting with the atom. 
The topology of these surfaces, determined by the coni- 
cal intersections, presents insight into the various atomic 
population and photon transfers from the maser field into 
the cavity field that can be produced by designing an ap- 
propriate path connecting the initial and the the chosen 
final states. Each path corresponds to a choice of the 
envelope of the pulses. In the adiabatic limit, when the 
pulses vary sufficiently slowly, the solution of the time- 
dependent dressed Schrodingcr equation follows the in- 
stantaneous dressed eigenvectors, following the path on 
the surface that is continuously connected to the initial 
state. We start with the dressed state | — ,0,0), i.e., the 
lower atomic state with zero photons in the cavity field. 
Its energy is shown in Fig. 2 as the starting point of 
the various paths. The paths shown in Fig. 2 describe 
accurately the dynamics if the time dependence of the 
envelopes is slow enough according to the Landau-Zencr 
m ill and Dykhne-Davis-Pcchukas [H ^ analysis. If 
two (uncoupled) eigenvalues cross, the adiabatic theorem 
of Born and Fock shows that the dynamics follows 
diabatically the crossing. This implies that the various 
dynamics shown in Fig. 2 are a combination of a global 
adiabatic passage around the conical intersections and lo- 
cal diabatic evolutions through (or in the neighborhood) 
of conical intersections of the eigencnergy surfaces [28l |. 

We consider the action of two smooth pulses, associ- 
ated with the Rabi frequencies G{t) and 0(i), which act 
on the two-level atom with a time delay t = d/v. Figure 
2 shows two examples of the adiabatic paths of different 
peak amplitudes of the Rabi frequencies and leading to 
two different photon transfers into the cavity field with- 
out changing the atomic population at the end of the 
interaction. Each of the two black paths (labelled (a) 
and (b)) corresponds to a sequence of two smooth pulses, 
shown in Fig. 3(a) and (c), of equal length Tint and dif- 
ferent peak Rabi frequencies fimax, Gmax , separated by 



+ 
to 

0) 

LU 




FIG. 2: First five eigenenergy surfaces (in units of 5) of as 
functions of G and SI for S — 2Am = —2Ac- The term G has 
been added to i?"** for clarity of display. The solid paths(a,b) 
correspond to adiabatic evolutions which start from | — ,0,0) 
state and end at | — , 1, —1) and | — , 2, —2) respectively. 



a delay such that the cavity pulse interacts before the 
maser pulse. 

For the path (a), the dynamics goes through the first 
intersection (on the f2 = plane) between the first and 
the second surfaces, but not the second intersection be- 
tween the second and the third surface. The crossing of 
the first intersection as G increases with $7 = 0, brings 
the dressed system into the second eigenenergy surface. 
Turning on and increasing the amplitude f2 (while G de- 
creases) moves the path across this surface. When the 
maser field decreases, the curve crosses another inter- 
section between the second and the third surface (with 
G = 0) that brings the system to the third surface, on 
which the path (a) stays until the end of the pulse fi. The 
transfer state is finally connected to the state | — , 1, — 1): 
there is no final transfer of atomic population, but one 
lum photon has been absorbed from the maser field and 
one ujc photon has been emitted into the cavity field at 
the end of the process. The path (b) allows the dynamics 
(on the r2 = plane) to go through the second intersec- 
tion, but not the third intersection. The next two inter- 
sections of the path (b) are located (on the plane G = 0) 
between the third and the fourth and between the fourth 
and the fifth eigenenergy surfaces and the system is fi- 
nally connected to the state | — , 2, —2): there is again no 
final transfer of atomic population, but two com photons 
have been absorbed from the maser field and two uJc pho- 
tons have been emitted into the cavity field at the end of 
the process. If the peak amplitudes are taken even larger 
such that n conical intersections (dynamical resonances) 
are crossed when G rises with f2 = and then n intcrsec- 
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Time (in units of 1/5) Time (in units of 1/5) 



FIG. 3: (a) and (c): Normalized Rabi frequencies with dif- 
ferent values of amplitudes for one-photon and two-photon 
transfers as a function of time, (b) and (d): Time evolution 
of the populations for one and two-photon transfers. 




1 2 3 1 2 3 



G/5 G/8 

FIG. 4: Contour plot of the populations P-,i,-i and P-,2,-2 
for one-photon and two-photon transfers as a function of nor- 
malized Rabi frequencies. 



tion are crossed when ft decreases with G = 0, the final 
state of the system will be \ — ,n, —n), i.e., the emission of 
n ujc photons into the cavity field and the absorbtion of 
n lom photons from the maser field, with no final atomic 
population transfer. The analysis of the eigenenergy sur- 
faces allows one to determine adapted amplitudes of the 
Rabi frequencies which will permit to transfer n photons 
into the cavity field in a robust way. 

IV. NUMERICAL SIMULATION 

The time evolution of the system for each family of the 
eigenenergy surfaces is given by the Schrodinger equation 

= H^'^imt) (22) 

The time dependence of the Rabi frequencies are delayed 
Gaussians of the form Figures 3(a), (c) show the pro- 
file of the Rabi frequencies as functions of time for one- 
photon and two-photon transfer with an interaction time 
(full width at half maximum) Tint = 66/(5 and a time 
delay t = 57/5. The adapted amplitudes of the Rabi 
frequencies for n-photon (n=l,2) transfer correspond to 
the paths (a) and (b) in Fig. 2. The condition for global 
adiabaticity \AM^c\Tint — = 33 » 1 is well satis- 

fied. Figures 3(b,d) present the time evolution of popula- 
tions calculated numerically by solving (|22|1 . The atom- 
maser-cavity system in the initial state | — , 0, 0) with the 
suitable forms of Rabi frequencies (Figs. 3(a,c)) evolves 
to the finale states | — , 1, —1) and | — , 2, —2) respectively 
with probabilities of P_ 1 _i = |(-, 1, -l|$(t/))|2 = 0.99 
andP_,2,-2 = |(-,2,-2|$(i/))|2 = 0.98. 



Fig. 4 displays the contour plot of the final popula- 
tion as a function of the normalized Rabi frequencies for 
one and two-photon transfers. The white regions repre- 
sent the adapted values of the Rabi frequencies for which 
the final probability of one and two-photon transfers are 
maximal. This figure shows that the bichromatic adia- 
batic passage is more robust with respect to the maser 
Rabi frequency than with respect to the cavity one. The 
reason comes back to the special structure of the dressed 
eigenenergy surfaces in Fig. 2. We can see that on the 
G = plane, between the first surface and the second 
one there isn't any intersection and between the ??th sur- 
face and its neighboring surfaces there are (n — 1) inter- 
sections, i.e. after the (n — l)th intersection there are 
no others. On the other hand on the = plane, as 
the value of G increases, the distance between neighbor- 
ing intersections decreases. In general, as the distance 
of conical intersections between neighboring surfaces de- 
creases, the robustness of the adiabatic passage is also 
decreased. 



V. DISCUSSIONS AND CONCLUSIONS 

Using the topological properties of dressed eigenenergy 
surfaces of the effective Hamiltonian of the atom-mascr- 
cavity system, we have determined adiabatic paths to 
transfer n photons from the maser field into the cavity 
field to generate a 71-photon Fock state. The realization 
of parameters satisfying the conditions of the proposed 
scheme appears feasible with progressive improvements 
to experiments with high-Q microwave cavities. In this 
analysis we have assumed that the interaction time be- 
tween the two-state atom and the fields is short compared 
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to the cavity lifetime Tcav and the atom's excited state 
lifetime Tatj i-G- ^int ^ Tcavi^at; which are essential for 
an experimental setup and avoiding decoherence effects. 
In the microwave domain, the radiative lifetime of circu- 
lar Rydberg states - of the order of Tat = 30ms - are 
much longer than those for noncircular Rydberg states. 
The typical value of the cavity lifetime is of the order 
of Tcav = 1 ms (corresponding to Q = 3 x lO'^) and the 
upper limit of interaction time is Tint = 100 jis (atom 
with a velocity of 100 m/s with the cavity mode waist of 
Wc = Q mm) m. 

The condition of global adiabaticity Gmax Tint ^ 1 for 
the typical value of Gmax ~ 0.15 MHz |23 is well satisfied 



(Gmax Tint ~ 15). Aualysis for additional conditions to 
go diabatically through conical intersections can be found 
in 
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